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Abstract 



We present a derivation of the gap equation for the crystalline color super- 
conducting phase of QCD which begins from a one-loop Schwinger-Dyson 
equation written using a Nambu-Gorkov propagator modified to describe the 
spatially varying condensate. Some aspects of previous variational calcula- 
tions become more straightforward when rephrased beginning from a dia- 
grammatic starting point. This derivation also provides a natural base from 
which to generalize the analysis to include quark masses, nontrivial crystal 
structures, gluon propagation at asymptotic densities, and nonzero temper- 
ature. In this paper, we analyze the effects of nonzero temperature on the 
crystalline color superconducting phase. 
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I. INTRODUCTION AND SUMMARY 



It is becoming widely accepted that at asymptotic densities the ground state of QCD 
with three massless quarks is the color-flavor locked (CFL) phase This phase features 

a condensate of Cooper pairs of quarks which includes ud, us and ds pairs. The CFL phase 
persists for finite masses, and even for unequal masses, so long as the differences are not too 
large It a ls° persists in the presence of a nonzero electron chemical potential /i e , so 

lone; not too large 0. In the absence of any interaction (and thus in the absence 

of CFL pairing) either a quark mass difference or a nonzero He pushes the Fermi momenta 
for different flavors apart, yielding different number densities for different flavors. In the 
CFL phase, however, the fact that the pairing energy is maximized when u, d and s number 
densities are equal enforces this equality [El. This means that if one imagines increasing 
either the strange quark mass m s or He, nothing happens until a first order phase transition, 
at which CFL pairing is disrupted, (some) quark number densities spring free under the 
accumulated tension, and a less symmetric state of quark matter is obtained . 

We can study much of the physics of interest by focussing just on pairing between massless 
up and down quarks with chemical potentials 

Hu = fi - 5/J, 

fi d = fx + Sfi. (1) 

For < 5ji < 5 Hi, the ground state is precisely that obtained for 5fi = In this state, 

red and green up and down quarks pair, yielding four quasiparticles with superconducting 



gap A pH|-[13|. And, the number density of red and green up quarks is the same as that of 
red and green down quarks. As 5/j, is increased from zero, this BCS state remains unchanged 
(and favored) because maintaining coincident Fermi surfaces maximizes the pairing and thus 
the gain in interaction energy. As 5fi is increased further, the BCS state remains the ground 
state of the system only as long as its negative interaction energy offsets the large positive 
free energy cost associated with forcing the Fermi seas to deviate from their normal state 
distributions. In the weak coupling limit, in which Aq/H <C 1, the BCS state persists while 
5h<5hi = A /V2 0,1- 

For 5fi ^> 5 Hi, only very weak pairing between like-flavor quarks is possible JR]. Near 



the unpairing transition, however, another phase intervenes. This is the "LOFF" state, 
first explored by Larkin and Ovchinnikov and Fulde and Ferrell in the context of electron 
superconductivity in the presence of magnetic impurities fll5| , |T6|| . Translating LOFF's results 
to the case of interest, the authors of Ref. || found that for S/j, ^ 8 Hi it is favorable to form 
a state in which the u and d Fermi momenta are given by Hu and Hd as in the absence 
of interactions, and are thus not equal, but pairing nevertheless occurs. Whereas in the 
BCS state, obtained for Sh < 6(/>i, pairing occurs between quarks with equal and opposite 
momenta, when <5/i > 5hi it is favorable to form a condensate of Cooper pairs with nonzero 
total momentum. This is favored because pairing quarks with momenta which are not 
equal and opposite gives rise to a region of phase space where each of the two quarks in a 
Cooper pair can be close to its Fermi surface, even when the up and down Fermi momenta 



1 



differ, and such pairs can be created at low cost in free energy.Q Condensates of this sort 
spontaneously break translational and rotational invariance, leading to gaps which vary 
periodically in a crystalline pattern. If in some shell within the quark matter core of a 
neutron star (or within a strange quark star) the quark chemical potentials are such that 
crystalline color superconductivity arises, as occurs for a wide range of reasonable parameter 
values, rotational vortices may be pinned in this shell, making it a locus for glitch formation 
0j. Rough estimates of the pinning force suggest that it is comparable to that for a rotational 
vortex pinned in the inner crust of a conventional neutron star, and thus may yield glitches 
of phenomenological interest ||. 

The authors of Ref. || studied crystalline color superconductivity in a simplified model 
with two flavors of quarks with chemical potentials ([!]) which interact via a four-fermion 
interaction with the quantum numbers of single gluon exchange. In the LOFF state, each 
Cooper pair has total momentum 2q with |q| « 1.25 fx. The direction of q is chosen sponta- 
neously. The LOFF phase is characterized by a gap parameter A and a diquark condensate, 
but not by an energy gap: the quasiparticle dispersion relations vary with the direction of 
the momentum, yielding gaps which vary from zero up to a maximum of A. The conden- 
sate is dominated by those regions in momentum space in which a quark pair with total 
momentum 2q has both members of the pair within ~ A of their respective Fermi surfaces. 
The gap equation which determines A was derived in Ref. [Q using variational methods, 
along the lines of Refs. [16|.[19| . This gap equation can then be used to show that crystalline 
color superconductivity is favored over no ud pairing for 5 fx < 5 1x2. Here, 5fX2 ~ 0.754A if 
the coupling is weak |T^,[T6|,[l9j and if there is no interaction in angular momentum J = 1 
channels M. For stronger coupling and for varying choices of interaction, 5fX2 changes |J. 

Crystalline color superconductivity is favored for 5fX\ < 5fx < 5fx 2 . As 5fx increases, one 
finds a first order phase transition from the ordinary BCS phase to the crystalline color 
superconducting phase at 5 fx = 5fi\ and then a second order phase transition at 5 fx = 5fx 2 
at which A decreases to zero. Analysis of the Ginzburg-Landau effective potential which 
describes physics near 5fX2 shows that A ~ (5 1x2 — 5fx) 1 ^ 2 for 5 fx — > 5 1x2 @. Because the 
condensation energy in the LOFF phase is much smaller than that of the BCS condensate at 
5fx = 0, the value of 5fX\ is almost identical to that at which the naive unpairing transition 
from the BCS state to the state with no pairing would occur if one ignored the possibility of 
a LOFF phase. For all practical purposes, the LOFF gap equation is not required in order 
to determine 5fx\. The LOFF gap equation is used to determine 5fi2 and the properties of 
the LOFF phase ||. For example, it determines the coefficients in the Ginzburg-Landau 
effective potential valid near 5fx 2 

The variational derivation of the gap equation for the crystalline color superconducting 
phase is somewhat cumbersome [01. One constructs a variational ansatz in which only 



x LOFF condenates have also recently been considered in two other contexts. In QCD with fi u < 0, 
fid > and fi u = —fid, one has equal Fermi momenta for u antiquarks and d quarks, BCS pairing 
between them, and consequently a (ud) condensate ]l7|]. If —u u and fid differ, and if the difference 



lies in the appropriate range, a LOFF phase with a spatially varying (ud) condensate results [17]. 
Suitably isospin asymmetric nuclear matter may also admit LOFF pairing, as discussed recently 
in Ref. El. 
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quarks within a "pairing region" are allowed to pair, minimizes the free energy with respect 
to all variational parameters (two per mode in momentum, color, flavor and spin space), and 
obtains a self-consistency relation which may then be solved to obtain A. The intricacy arises 
from the fact that the definition of the boundary of the pairing region involves A itself. A 
derivation in which one simply makes an ansatz for the quantum numbers of the condensate 
and then "turns a field-theoretical crank" and sees this intricate result emerge would be 
helpful both by virtue of being more straightforward and because the use of variational 
methods to obtain a gap equation is by now less familiar to many readers. We provide such 
a diagrammatic derivation here. 

Furthermore, and as we explain at appropriate points in our presentation of the derivation 
in Sections II, III and IV, many generalizations are amenable to analysis using the formalism 
we present here: 

• In Section V, we include the effects of nonzero temperature. We calculate the critical 
temperature T c above which the crystalline color superconducting condensate vanishes 
and show that for 5jx — > 5fi2, T c — > 0.39 A, as previously known |L9|. 



It should also be straightforward to generalize the analysis to include three flavors of 
quarks with differing masses, and thus to study the crystalline color superconducting 
phases expected where either fi e or m s is just larger than that at which the CFL phase 
is lost ||. This, not the toy model we analyze here, is the case of physical interest. 

As in Ref. 0, we restrict our attention here to the simplest possible "crystal" struc- 
ture, namely that in which the condensate varies like a plane wave. Wherever this 
condensate is favored over the homogeneous BCS condensate and over the state with 
no pairing at all (i.e. where Sfx\ < 5li < 8^2) we expect that the true ground state of 
the system is a condensate which varies in space with some more complicated spatial 
dependence. The formalism we set up can be generalized to derive gap equations for, 
and hence to analyze and compare, condensates with arbitrary crystal structures in 
order to learn which one is favored. 

The diagrammatic analysis we present uses a point-like interaction between quarks, 
but the formalism is easily generalized to treat the exchange of a propagating gluon, 
as appropriate at asymptotically high densities. Even if such analyses, pioneered at 
5 ix = in Ref. [^(| and since studied in considerable detail by many authors 0, are 
to date of quantitative value only at inaccessibly high densities pi| , it would be very 
interesting to see crystallization occurring in a controlled analysis beginning directly 
from the QCD Lagrangian. 



II. THE GAP EQUATION FOR CRYSTALLINE COLOR 
SUPERCONDUCTIVITY 

In the ordinary BCS phase, pairing between quarks with momentum p and — p is de- 
scribed in the standard Nambu-Gorkov formalism by introducing an eight-component field 
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such that, in this basis, the inverse quark propagator takes the form 



s-\p) 



^ + /i7o A(p) 
A(p) (j> - /i 7o ) T 



(3) 



Here, A = 7 A^7 and A is a matrix with color, flavor and Dirac indices which have all been 
suppressed. The diagonal blocks correspond to ordinary propagation and the off-diagonal 
blocks reflect the possibility of "anomalous propagation" in the presence of a diquark con- 
densate (?/>(x)-?/>(x)) oc A. 

In the crystalline color superconducting phase [|T|,|l6|,|] , the condensate is made up of 



pairs of u and d quarks with momenta such that the total momentum of each Cooper pair is 
given by 2q with |q| « 1.25 fx. The direction of q is chosen spontaneously. Such a condensate 
varies periodically in space, with wavelength 7i7|q|: 



(^(x)^(x)) oc Ae 2iqx 



(4) 



The spatial dependence (|]) is only the simplest possible choice. Wherever is favored 
over both the ordinary BCS state and the state with no pairing at all, we expect that the 
true ground state of the system will include Cooper pairs with their respective |q|'s taking 
on the same, energetically favored, value, but choosing one of several spontaneously selected 
directions. The result would be a condensate which varies in space like a sum of plane 
waves. For example, a cubic crystal arises as a sum of six plane waves. The favored crystal 
structure for the crystalline color superconductor is not known. In this paper, as in Ref . , 
we only consider the simplest possibility (|). 

Although we expect that crystalline color superconductivity occurs whenever the mass 
difference or chemical potential difference between any two flavors of quarks is just larger 
than the maximum value which the standard BCS state can tolerate, for concreteness we 
shall follow Ref. |J and only consider pairing between massless u and d quarks with chemical 
potentials ([]]). In the condensate (f|), u quarks with momentum p + q pair with d quarks 
with momentum — p + q. To describe this, we must use a modified Nambu-Gorkov spinor 
defined as 



*(p,g) 



td(p - q) 
tJ(-p + q) 
V ^u(-p - q) J 



(5) 



Note that flavor indices are now explicit, which will be convenient below. The central change 
we have made in going from (0) to (0) is to modify the momentum dependence. Note that 
by q we mean the four- vector (0, q). The Cooper pairs have nonzero total momentum, 
but the ground state condensate (f|) is static. The change from (fj) to (||) can be seen as a 
change of basis. In the presence of a crystalline color superconducting condensate, anomalous 
propagation does not only mean picking up or losing two quarks from the condensate. It 
also means picking up or losing momentum 2q. If we tried to describe this using the original 
basis (0), the inverse quark propagator would no longer be diagonal in momentum space. 
The new basis (|J) has been chosen so that the inverse quark propagator in the crystalline 
color superconducting phase is diagonal in p-space and is given by 
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S-\p,q) 





-A(p, -q) 






- i + Vdio 
o 

A(p,g) 



-A(p, -q) 
A(p,g) 
(p'-fi-Vdlof 

# + i - Vulof 



(6) 



2p is the relative momentum of the quarks in a given pair, and is different for different pairs. 
In the gap equation below, we shall integrate over p and p, as we sum the contribution of all 
pairs. 2q is the center of mass momentum of every pair in the condensate; it is a constant 
and thus will not be integrated over. It is convenient to denote flavor indices explicitly 
in (||) because we are describing the situation where \i u ^ fid- It is straightforward to 
introduce different quark masses in Eq. (^j), but then the calculations become more involved 
and we therefore defer this to a future publication. Note that the condensate is explicitly 
antisymmetric in flavor. Color and Dirac indices remain suppressed. As desired, the off- 
diagonal blocks describe anomalous propagation in the presence of a condensate of diquarks 
with momentum 2q. The choice of basis we have made is analogous to that introduced 
previously in the analysis of a crystalline quark-antiquark condensate [22]. This work also 
points the way toward the generalization of ([|) needed to handle a condensate which varies 
in space like a cubic crystal rather than the plane wave (||). 

We wish to obtain the gap by solving the one-loop Schwinger- Dyson equation, given by 



S-\k, q) - S \k, q) = ig 2 J t^j^Cp, D&(k - p) 



(2tt)< 



(7) 



Here, D^ B is the gluon propagator, S is the full quark propagator, whose inverse is given by 
(§), and Sq is the fermion propagator in the absence of interaction, given by S with A = 0. 
So looks unusual, because it depends on both k and the "offset" q. This is a consequence of 
our choice of basis (H), and would be a legitimate if perverse way to describe noninteracting 
fermions. This choice of basis is natural in the crystalline color superconducting phase. The 
vertices are defined as follows: 



/7„A A /2 
l,\ A /2 


V o o 






-(7^/2) G 










\ 



(8) 



-( 7 ^/2f y 



Note that in Eq. (0) we have chosen to work in Minkowski space. We will continue to work 
in Minkowski space until we obtain the gap equation itself, which we will then write in 
Euclidean space for computational convenience. 

We defer the analysis of the crystalline color superconducting phase at asymptotically 
high densities to future work. In this paper, as in Ref. 0, we choose to caricature the inter- 
action between quarks as a point-like four-fermion interaction with the quantum numbers 
of single-gluon exchange. This means that in ([/]), we make the replacement 



9 2 D 



AB 



(9) 



where G, normalized as in Ref. 0, is a dimensionful coupling constant which parametrizes 
the strength of the interaction between quarks. Reasonable choices for G, motivated by zero 



5 



density hadron phenomenology, yield a BCS gap on the order of 100 MeV at fl = 400 MeV 
in the absence of any chemical potential difference 8 fx 0. 

Once we have removed the gluon propagator from ([?[), we see that the right-hand side is 
independent of k. The left-hand side must therefore be independent of k as well, meaning 
that A in (^) must be independent of p. A does depend on the common momentum of all the 
Cooper pairs, 2q. Choosing an ansatz for A is straightforward once we have understood that 
it must be independent of p. Single gluon exchange is attractive in the color-antisymmetric 
(3), flavor antisymmetric, Lorentz scalar or pseudoscalar channels. Instanton effects favor 
the scalar condensate, and we therefore make the ansatz 

A a P{p,q) = e al33 C-f 5 A (10) 

for the gap matrix, where a and j3 are color indices, running from 1 to 3, and C = 27072- 
A has no remaining indices and all the matrix structure has now been written explicitly. 
A does depend on |q|, although we do not denote this dependence explicitly, but does not 
depend on p or the direction of q. 

After some algebra (essentially the determination of S given S 1-1 specified above), and 
upon suitable projection, the Schwinger- Dyson equation (|7|) reduces to a gap equation for 
the gap parameter A given (in Euclidean space) by 

, =2C [^P_ 1^ (n) 

J (2vr)^2_ 4 [ ( | p |2_ ( ^ o + W)(/i2 _| q |2 ) + (p . q + /i( ^ o + 5/i)) 2 ] ^ ) 

where w = |p| 2 — |q| 2 — (ip + 8fx) 2 + fx 2 + A 2 . 

We analyze the gap equation (|TJ) in Section IV. It will turn out to be close to, but 
not identical to, that derived in Ref. ||. The difference is that here we have kept the 
contributions of particles, holes, and antiparticles in the gap equation, whereas in Ref. 
@ only particle-particle and hole-hole pairing was considered. Pairing in the crystalline 
color superconducting phase is dominated by those pairs in which both particles or both 
holes in a pair are near their respective Fermi surfaces. Indeed, it is the fact that such 
pairs exist even at nonzero 8 fx as long as |q| > 8 fx which explains why the crystalline 
color superconducting phase may be favored in the first place. We therefore expect that 
neglecting the contributions of the antiparticles in the gap equation, as was done in Ref. || , 
should be a good approximation. Demonstrating this requires complicating the gap equation 
considerably at first, although it does eventually simplify as will be shown by the end of 
Section III. We shall see in Section IV that once the contributions of antiparticles have been 
eliminated, the gap equation we derive here agrees with that of Ref. || . 



III. ELIMINATING ANTIPARTICLES 

In order to eliminate the (small) contribution of the antiparticles on the right-hand side 
of the gap equation, we shall need the projectors p3| , p4 



P+(P) = l±±l 

P-(P) = ^ (12) 



6 



with a = 707, where P+ projects onto particle states and P_ projects onto antiparticle states. 
This allows us to separate the ansatz for A into those parts which include antiparticle pairing 
and those which do not. We replace the ansatz ( |T0| ) by 

A a ?(p, -q) = e a/33 C 7 5 [AiP+ip - q)P + {p + q) + A 2 P_(p - q)P_{p + q) 

+ A 3 P + (p - q)P_(p + q) + A 4 P_(p - q)P + {p + q)\ . (13) 

Here, A]^ 3 4 are four (potentially different) gap parameters whose meaning we now explain. 
To understand each of the terms in fljjD , note that, for example, 

C 7 5 P+(p - q)P+(p + q) = P£(- P + q)C 7 5 P + (p + q) . 

Thus, Ai describes pairing between particles (and not antiparticles) with momenta p+q and 
— p + q. Similarly, A 2 describes antiparticle-antiparticle pairing, and A 3 and A 4 describe 
particle-antiparticle pairing, which is only possible for q 7^ 0. (For q = 0, the only projectors 
which occur are P+{p) and P-(p), and P + (p)P_(p) = 0.) 

With our point-like interaction, A must be independent of p. This requires 

Ai = A 2 = A 3 = A 4 = A , (14) 



which restores the simple ansatz (|10J) . It may seem perverse, but we now derive coupled 
gap equations for Ai^.3,4, without assuming that they are equal. We do so for two reasons. 
First, in a future publication, we plan to restore the gluon propagator. In this context, A 
is not independent of p and Ai 2,3 4 therefore need not all be the same. The exercise below 
therefore lays the groundwork for this future calculation. Second, and in the present context, 
we wish to eliminate all terms on the right-hand side of the gap equation which depend on 
A2,3,4, as they make only a small contribution. The reader not interested in details of this 
derivation can safely skip to Eqs. ( pi| ) an d ( |2"2"D and the discussion that follows them. 
The Schwinger-Dyson equation (]?]) (using Eq. (§)) is an equation for the matrix A: 

A(k,q) = (l^ T S A2 {p,q) (yy) (15) 

where S42 (k, q) is the (4,2)-component of the fermion propagator found by inverting Eq. (^j). 
After some algebra, we find 



+ W° ~ A(p, q)(i> + i- /i u 7°) A(p, q) 



-1 



(16) 

Upon inserting the ansatz ( ]T3|) for A(p, q), we can rewrite (|l^) in terms of the gap parameters 
Ai.,2,3,4. In order to display the resulting expressions, we must first define: 

, . Aisin 2 (/5/2) + I/2A3COS/3 . A 4 cos 2 (/3/2) - 1/2 A 2 cos/3 

A = p + fi d -\p-q\- Ai — — ■ — : A 4 ; — ■ — , 

Po ~ + |p + q| Po - - |p + q| 

A 2 sin 2 (/3/2) + l/2A 4 cos/3 A 3 cos 2 (/3/2) - l/2Ai cos/3 
i i A3 ; ; (17 ) 

Po - v-u - |p + q| Po - + |p + q| 

_ A 2 A 4 \ 

po- ^u- |p + q|y 
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Here, the angle (3 is defined as the angle between the up quark momentum q + p and the 
down quark momentum q — p and is therefore given by 

cos/3= (q + p) • (q^p) • (18) 
With these definitions, S i2 (p,q) becomes 

S 4 2(p,q) = e^ 3 C l5 T(p,q) (19) 

with 

T(p, q) = ~ \xs=B^=m^9 ( po -^ + | p + q | P -(P + 9) P "(P - 9) 

+ 4 1 , \ p +(p+q)P-(p-q)) 

+ ^W-B)Cco S/3 ( po _ /iu + 3 |p + q| ^^ + q)P + iP - q) 

+ A '| ~ , P+(j> + q)P+(p - q)) (20) 

A| A " \P+{p + q)P+{p - q)P-{p + q) 



Pq- Vu - |P + q| 



Aisin 2 f + A 3 cos 2 f A 2 sin 2 f + A 4 cos 2 f , 



■P- (p + g) + — 2 - — r f-P+ (p + g) 



Noting that 



Po-Ai« + |P + q| po-/^«-|p + q 



3 



we obtain the gap equation 

AxP + (A; + q)P + (k -q) + A 2 P_(/c + q)P_{k - q) 
+A 3 P + (k + q)P_(k -q) + A A P_(k + q)P + (k - q)= -2iG J 7^7^ q)l" ■ (21) 



Upon setting all the A's equal as in flUD, the gap equation ( pi]) yields 

A = -2 ? G/^TrT(p,g). (22) 
When written explicitly, this is Eq. (ITT]). 



Using the definition of T(p,q), we can easily identify the contributions of the four A's 
to the right-hand side of the gap equation. Since we expect that the Ai terms, which 
describe particle-particle and hole-hole pairing, will give the dominant contributions to the 
gap integral, we now eliminate all terms in T(p, q) which depend on A 2i3j4 . This means 
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that the contribution of the antiparticles to the right-hand side of the gap equation has been 
eliminated. Note that we are not setting A2.34 = 0- With a point-like interaction, all the A's 
are in fact equal as in (0), and the left-hand side of the gap equation (|2T| ) is /^-independent 
and equal to A. The point is that the contributions of those terms on the right-hand side 
of (|2T| ) in which A2.34 appear must be small, and we can therefore neglect them. In other 
words, once all integrations have been completed on the right-hand side, we would find that 
^2,3,4 only occur multiplied by quantities which are small if A//i is small. Dropping A2,3,4 
on the right-hand side before integration (but keeping them on the left-hand side) should 
therefore be a good approximation. The resulting gap equation can be written (in Euclidean 
space) as 

r f d *P 2 A sin2 f 

where £?i )2 (p) are defined as in Ref. |§: 



£a(p) = + Sfx + \ (|p + q| - |p - q|) + \ J(|p + q| + |p - q| - 2/i) 2 + 4A 2 sin 2 f 



E 2 (p) = - 5fi - \ (|p + q| - |p - q|) + I^Ip + q| + |p - q| - 2/x) 2 + 4A 2 sin 2 f (24) 

and f3 is defined in Eq. (|Tj|). As we describe below, we have confirmed explicitly that the 
gap equation (p^) is a good approximation to Eq. QTT|). 



IV. ANALYZING THE GAP EQUATION 



The energies E\ and E 2 given in Q2"^] ) arise in Ref. || . There, we deduced that the right- 
hand side of the gap equation must be taken to vanish in those regions of p-space where 
either £a(p) or E 2 (p) is negative via the following argument. In the region where £a(p) < 0, 
it is free-energetically favorable to have unpaired w-quarks rather than pairs. Similarly, in 
the region where ^(p) < 0, it is free-energetically favorable to have unpaired <i-quarks 
rather than pairs. Because quarks do not pair in these "blocking regions" of momentum 
space, these regions do not contribute to the gap equation, which becomes an integral over 
those regions of momentum space wherein pairing occurs. If A is set to zero, the blocking 
regions are simply described. They are the regions in p-space where the w-quark state with 
momentum p + q is within the u Fermi sea while the d-quark state with momentum — p + q 
is outside the d Fermi sea, or vice versa. In the presence of a nonzero A, the boundaries of 
the blocking regions are given by £?i(p) = and ^(p) — and therefore depend on A and 
are not simply determined by the locations of the noninteracting Fermi surfaces. The result 
of this analysis, presented in Ref. ||, is a variational procedure in which the boundaries of 
the blocking regions, and thus the specification of the variational ansatz itself, depend on 
the gap A, which is in turn obtained by solving a gap equation whose integrand is restricted 
by hand to vanish within said blocking regions. 

In contrast to the intricacy of the variational approach, the physics of the blocking regions 
emerges from a completely straightforward analysis of the gap equation in the form we have 
derived above, namely Eq. fl23|). We simply do the po integral by contour integration. There 
are two poles, both of which lie on the imaginary axis. Let us close the contour in the upper 
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half plane. If E\ > and E 2 > 0, we pick up the pole at po = %E\ which has a residue 
proportional to + E 2 ). If E\ < and -E2 > 0, both poles are in the lower half plane, 

and the right-hand side of the gap equation vanishes. If E\ > and E 2 < 0, both poles are 
in the upper half plane, the residues from the two poles cancel, and the right-hand side of 
the gap equation again vanishes. (If we close the contour in the lower half plane, we obtain 
the same result upon noticing that we encircle no poles if E\ < and E 2 > and two poles 
with cancelling residues for E\ > and E 2 < 0.) Thus, upon doing the po integration we 
obtain the gap equation of Ref. 0: 



2G 



d 3 p 2 sin 2 § 
P ev (27r)3 J B 1 (p) + E 2 (p) 



2G[ J% - ^ (25) 

( 27r ) M\p + q| + |p - q| - 2/i) 2 + 4A 2 sin 2 § 



where the "pairing region" V in p-space is given by 

V = {p I Ei(p) > and E 2 (p) > 0} . (26) 

Thus, a trivial exercise in residue calculus has reproduced the blocking regions, excluding 
from the gap equation those regions in momentum space where Ei(p) or ^(p) is negative. 
Note that because -E'i(p) + -^(p) — 0, as can be seen from the definitions (0), there is 
no value of p for which both E\ and E 2 are negative. Note also that the gap equation is 
dominated by those regions in momentum space where -Ei(p) + E 2 (p) is as small as possible, 
where the integrand in fl2"5| ) is of order 1/A. These values of p are such that both members 
of a LOFF pair have momenta close to (within ~ A of) their respective Fermi surfaces. 
That is, |p + q| is within A of fi u and |— p + q| is within A of [id- 

For completeness, we sketch the analysis of (|TT|), in which the contributions of antiparticle 
pairing to the gap equation have not been eliminated. The denominator of the integrand in 
(|lTJ) is a fourth order polynomial in p , so the gap equation can be rewritten as 

A = 2G \ — v — 7 ^ ^ . 27 

J (2vr)4 (p _ , Pl(p) ) ( po + iPa ( p )) fa + iPl(p )) fa _ ,p 2 ( p )) 

The analytical expressions for the poles -Pi(p), -^(p), A(p), and -P2(p) are complicated and 
uninformative. However, we have checked that for reasonable choices of parameters, the 
numerical values of -Pi(p) and ^(p) are very close to -Ei(p) and E 2 (p) and those of -Pi(p) 
and P 2 (p) are very close to the antiparticle energies 



^i(p) = -fy + i(|p + q|-|p-q|) + iV(|p + q| + |p-q| + 2/i) 2 + 4A 2 sin 2 ^ 



E 2 (p) = + 8fi - \ (|p + q| - |p - q|) + \sj(\p + q| + |p - q| + 2ft) 2 + 4A 2 sin 2 f . (28) 

The analysis of ( |2"T| ) is analogous to that of (^). Wherever both P 1 and P 2 are positive, 
the pole at p = iPi contributes, with residue proportional to l/(Pi + P 2 ). The integral 
is dominated by the region where Pi + P 2 is close to zero. There is also a contribution 
from the pole at po = iP 2 , but the residue of this pole is nowhere large. For this reason, 
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and because Pi and P2 are numerically very close to E\ and E2, we find that (0), and 
thus (^) which was derived variationally in Ref. [§, is a good approximation to ([27]). For 
/i = 400 MeV and G chosen such that the BCS gap at = is 100 MeV, we take Sfi within 
the range where the crystalline color superconducting phase is favored || and choose a value 



of A which solves fl23|). We then find that the right-hand sides of fl23|) and ( |2T|) differ by 
about 20%. The discrepancy vanishes in the weak-coupling limit. (Note that eliminating 
the contribution of the pole at po = 1P2 does not by itself reduce (^) to (|2"B|), although it 
does change the 20% discrepancy to a 1% discrepancy. Eliminating all the contributions of 
the antiparticles is more subtle, as we have seen.) 

Although the expressions for the P(p)'s are much more complicated than those for the 
i?(p)'s, the gap equation ( |27|) actually turns out to be more easily solvable (by Mathematica) 
than (p3[), as we now explain. The blocking regions in (|27| ) are regions wherein either P\ or 
P2 is negative, and are therefore bounded by surfaces on which Pi or P 2 vanishes. Within 
these regions, there is no contribution from the Pi and P2 poles. The simplification which 
occurs in fl27|) is in the explicit expressions for the boundaries of the blocking regions. We 
simply set po = in the denominator of Eq. ( |TT| ) and for each value of |p| we solve for cos 9, 
the angle between p and q: 



2flSn±yJ(\p 


2 _ 


q 


2 — 5jJ? + jJL 


2 - 


-A 2 ) 2 -4(| P 


2 -<5/i 2 )(/i 2 - 


q 


2 ) 


2 


pl 


|q 





cos9 = — — 5L^J LZ! : 1 11 1 — LZL - i ■ ( 29 ) 



This allows one to implement the fact that the p-integral is to be taken over all of p-space 
except for the blocking regions via explicitly specified limits on the cos 9 integral. In contrast, 
even though Ei and E2 are simpler than Pi and P2, the dependence of the denominator in 
(p3|) on cos 9 is more complicated than that in (|TT|), and the blocking regions can only be 
specified explicitly as roots of a quartic polynomial. 

The calculation of 5fi2 of Ref. || follows directly from Eq. (|23|), as does the value of 
|q|, which is given to a very good approximation by that for 5fi — > 5fi2- (One finds 8^2 
by seeking the largest value of <5/i for which there is a choice of |q| which yields a nonzero 
solution A to the gap equation (^3|) 0.) Our rederivation of (^3|) has several merits. First, 
as it begins with a Schwinger-Dyson equation rather than a variational wave function, it 
may appear more familiar. Second, the emergence of blocking regions is straightforward. 
Third, it is the basis for many generalizations: Quark masses can easily be introduced in (^j), 
and the analogue of ( pTTj ) can then be derived. The gluon propagator need not be replaced 
by a point-like interaction. This opens the way to a treatment of color superconductivity at 
asymptotically high density. Also, nontrivial crystal structures can be analyzed beginning 
with a Nambu-Gorkov propagator which admits "anomalous propagation" in which 2q of 
momentum is gained or lost, for several different values of the vector 2q. Finally, the 
generalization to nonzero temperature is straightforward. To this we now turn. 



V. CRYSTALLINE COLOR SUPERCONDUCTIVITY AT NONZERO 

TEMPERATURE 

We can now derive the gap equation for the LOFF state at nonzero temperature. We 
begin with Eq. (p3|). Using the standard formalism, we obtain the nonzero temperature gap 
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equation by converting the po integral into a sum over Matsubara frequencies. That is, with 
uj n = (2n + l)vrT, using the prescription 



Po - 

we obtain the following equation: 

d 3 p 



A = 2G 



-2 A sin" 



uj„ and 



'i ) * t 



dpo 
2tt 



1 



(2vr)3 V 2 / nf^oK-^i(p))K + ^ 2 (p)) 

The sum may be evaluated by converting it into a contour integral: 



(30) 



E 



[UJn 



dz 



iE 1 )(u n + iE 2 ) 2m Jc T(e z / T + 1) {z + E{){z - E 2 ) 



where the contour C encircles the imaginary axis.Q We may deform the contour so that it 
encircles, with negative orientation, the poles off the imaginary axis. This gives 



dz 



+ 



2m Jc T(e*/ T + 1) (z + E 1 )(z - E 2 ) T(e E */ T + 1) E 2 + E x T(e~ E ^ + 1) E 2 + E x 



1 



tanhl^ 
2T 



2T{E X + E 2 

Upon using these identities in the finite-temperature gap equation Eq 



tanh ( —^i 
2T 



Of), we obtain: 



2G 



d 3 p 



2sin 2 (f) 



(27r) 3 E 1 (p)+E 2 (p) 



tanh 



' £i(p) 
2T 



+ tanh 



' e 2 (pY 

2T 



(31) 



Note that here the integration is performed over all of p-space: there are no blocking regions 
at nonzero temperature. The blocking regions emerge in the limit T — >• as follows: if 
£i(p) > and E 2 (p) < 0, then tanh(^i^) 1 while tanh(^2i) -> -1, and the integrand 
vanishes. The same result holds if Ei(p) < and ^(p) > 0. If, however, -Ei(p) > and 
E 2 (p) > 0, 



tanh 



EM 
2T 



+ tanh 



E 2 (p) 
2T 



reproducing the zero temperature gap equation Eq. (P5|). Note that even if our goal were 
just to understand physics at T = it may be of practical value to do calculations at several 
nonzero values of the temperature and then extrapolate to T = 0. The reason is that at 
T = 0, specifying the boundaries of V of (p6|) , and thus the limits of integration, can be a 
numerical challenge. At any nonzero temperature, instead, no limits of integration need be 
specified. The tanh factors impose the required limits as T gets small. 

We shall present T / results for parameters chosen as in Fig. 4 of Ref . , which we 
first recapitulate. We specify the four-fermion interaction by choosing the cutoff parameter, 



2 Note that 



T{e^ T +l) 



has simple poles with unit residue at z 
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FIG. 1. The gap A as a function of temperature T, at <5/i = 5 fa. At zero temperature, 
A = 7.8 MeV= 0.195A . The gap vanishes above T c = 3.42 MeV. 



defined in Ref. @, to be A = 1 GeV and requiring G to be such that the BCS gap is A = 
40 MeV at 5fi = 0. We choose fl = 400 MeV, and explore different values of 5fi. At T = 
||, we find nonzero solutions to the gap equation (^) for 5fi < 6 fa = 0.744A . Above 8 fa, 
no pairing between u and d quarks is possible. The crystalline color superconductor phase 
has lower free energy than the ordinary BCS phase as long as Sfi > 5 fa = 0.710A , where 
a first order phase transition occurs. A precise determination of 5 fa requires expressions 
for the free energy of both phases. The free energy of the crystalline color superconductor 
phase could be obtained from the gap equation along the lines described in Section 4.3 of 
Ref. ||. However, 5 fa is well approximated by the Sfi at which the BCS and unpaired states 
have equal free energy, which turns out to be 0.711A ||. At 5fi = 8 fa, A = 7.8 MeV 0. 
For 5/i — > 5 fa from below, A vanishes like (5 fa — Sfi) 1 ^ 2 ||. The window 5 fa < 5fi < 5 fa 
widens if the interaction includes attraction in the spin-one channel. We expect this window 
to widen at asymptotic density, where quarks interact by exchanging a propagating gluon. 

In Fig. 1, we show the dependence of A on the temperature T at 5fi = 5 fa. We find that 
the critical temperature above which the crystalline color superconductivity is lost is T c = 
3.42 MeV, corresponding to T c = 0.44A(T = 0). In Fig. 2, we plot both T c and A(T = 0) 
as functions of 5 fa for 5 fa < 5fi < 5 fa. We find that the ratio T C /A(T = 0) changes little, 
decreasing from 0.44 at 5 fa to 0.39 for 5fi — > 5 fa. This agrees with the previously known 
result that T C /A(T = 0) — > y / 3/27r 2 for 5fi — > 5 fa Hl9fl . There are two a priori reasons why 
one may have questioned whether T c in the crystalline color superconducting phase would 
turn out to be proportional to A(T = 0). First, this phase is in fact gapless. There are 
directions in momentum space (which intersect the boundaries of the blocking regions) for 
which gapless excitations exist at zero temperature. One may wonder whether the presence 
of these gapless modes, which can be excited at arbitrarily low temperature, might lower T c . 
Second, the condensation energy in the crystalline color superconductor phase is of order 
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the ordinary BCS phase is favored. 

/i 2 A 4 /Aq |§, whereas that in the ordinary BCS phase is of order /2 2 Aq. One may therefore 
wonder whether the T c for crystalline color superconductivity scales differently with A. It 
turns out, however, that the simplest expectation holds true: A is the gap in the fermion 
spectrum in directions in momentum space along which pairing is maximized and destroying 
the condensate therefore requires a temperature T c which is of order A. 
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